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We consider a pair of parallel straight quantum waveguides coupled lat- 
erally through a window of a width i in the common boundary. We 
show that such a system has at least one bound state for any t > . 
I/-) ■ We find the corresponding eigenvalues and eigenfunctions numerically 

using the mode-matching method, and discuss their behavior in sev- 
eral situations. We also discuss the scattering problem in this setup, 
in particular, the turbulent behavior of the probability flow associated 
with resonances. The level and phase-shift spacing statistics shows that 
in distinction to closed pseudo-integrable billiards, the present system 
is essentially non-chaotic. Finally, we illustrate time evolution of wave 
packets in the present model. 

X 



I Introduction 

Spectral and scattering properties of quantum particles whose motion is confined 
to nontrivial subsets of M n represented until recently rather textbook examples 
or technical tools used in proofs. There are several reasons why these problems 
attracted a wave of interest in last few years . The most mathematical among them 
stems from the observation that, roughly speaking, one can choose the region in 
such a way that the spectrum of the corresponding Neumann Laplacian coincides 
with a chosen set 1,2 ; of course, the boundary of such a region may be in general 
rather complicated. 

On the other hand, even regions with nice boundaries may exhibit various unex- 
pected properties manifested, for instance, in spectra of the corresponding Dirichlet 
Laplacians. A prominent example is the existence of bound states, i.e., localized 
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solutions to the free Schrodinger equation, in infinitely stretched regions such as 
bent, branched or crossed tubes of a constant cross section — see, e.g., Refs. 3-7; 
more references are given in the review paper 8 . 

1.1 Quantum wire systems 

A strong motivation to study such bound states and related resonance effects 9-11 
comes from recent developments in semiconductor physics, because they can be 
used as models of electron motion in so-called quantum wires, i.e., tiny strips of 
a very pure semiconductor material, and similar structures. Let us briefly recall 
key features of such systems; more details and a guide to physical literature can be 
found in Ref. 8. 

Characteristic properties of the semiconductor microstructures under consider- 
ation are small size, typically from tens to hundreds of nm , high purity, which 
means that the electron mean free path can be a few /im or even larger, and crys- 
tallic structure. In addition, boundaries of the microstrucures consist usually of an 
interface between two different semiconductor materials; the electron wavefuction 
are known to be suppressed there. 

Behavior of an electron in such a "mesoscopic" system structure is, of course, 
governed by the many-body Schrodinger equation describing its interaction with the 
lattice atoms including the boundary, external fields, and possible impurities. The 
mentioned properties allow us, however, to adopt several simplifying assumptions. 
As we have said the mean free path is typically two or three orders of magnitude 
greater than the size of the structure; hence the electron motion can be assumed in 
a reasonable approximation as ballistic, i.e., undisturbed by impurity scattering. 

The most important simplification comes from the crystallic structure. The one- 
electron Hamiltonian as a Schrodinger operator with a periodic potential exhibits 
an absolutely continuous spectrum — see Ref. 12, Sec. XIII. 16 — in the solid-state 
physics language one says that the electron moves in the lattice as free with some 
effective mass m*. The latter changes, of course, along the spectrum but one can 
regard it as a constant when we restrict our attention to the physically interesting 
part of the valence band; recall that its value may differ substantially from the true 
electron mass, for instance, one has m* = 0.067 m e for GaAs which is the most 
common semiconductor material used in mesoscopic devices. 

This property together with the wavefunction suppression at the interfaces 
makes natural to model electrons in a quantum wire system as free (spinless) par- 
ticles living in the corresponding spatial region with the Dirichlet condition on its 
boundary; an interaction term must be added only if the whole structure is placed 
into an external field. This is the framework in which the mentioned curvature- 
induced bound states and resonances were studied. However, the physical conclu- 
sions one can draw from it are not restricted to mesoscopic devices: the results are 
useful for description of other "new" quantum systems 13 and provide fresh insights 
into the classical theory of electromagnetic waveguides 14-16 . 
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1.2 Motivation of the present work 

Apart from these practical reasons, the curvature-induced bound states provide at 
the same time a warning example showing that an intuition based on semiclassical 
concepts may fail when dealing with quantum systems. It is well known, for instance, 
that low-dimensional Schrodinger operators have bound states for an arbitrarily 
small coupling constant as long as the potential is not repulsive in the mean and 
decays sufficiently fast at infinity 17 ' 18 . A common wisdom, however, is that this is 
rather an exception, and that the number of bound states of a quantum system is 
at least roughly proportional to the classically allowed volume of the phase-space. 
The waveguide systems in question illustrate that this is not true, because they can 
exhibit in principle any number of bound states, while having no closed classical 
trajectories with the exception of an obvious zero-measure set. 

In this paper we are going to consider another example of that kind, this time 
consisting of two straight parallel quantum waveguides. We suppose that they have 
a common boundary which has a window of a width I allowing the particle to leak 
from one duct to the other. This an idealized setup for several recently studied 
quantum-wire systems — see Refs. 19-22. Using a variational argument we shall 
show that such a system has always at least one bound state, i.e., an isolated 
eigenvalue below the threshold of the continuous spectrum. Moreover, the system 
can have any prescribed number of bound states provided the window width is chosen 
large enough. These conclusions follow from simple estimates; however, they tell us 
nothing about the corresponding wavefunctions and more detailed dependence of the 
bound-state energies on the parameters. To this aim, we shall formulate in Section 4 
a method to solve the problem numerically using the mode-matching technique. In 
particular, we shall discuss how the first eigenvalue emerges from the continuum as 
the window opens. 

Interesting properties of the system are not exhausted by this. The coupling be- 
tween the wavefunctions in the "arms" and the connecting region allows the particle 
to tunnel betwen different transverse modes, so the scattering matrix is nontrivial; 
one naturally expects it to have resonances with properties similar to those of the 
bound states below the first transverse-mode energy. It is not easy to prove the 
existence of such resonances, because there is no natural parameter in the problem 
which would make it possible to tune the intermode coupling. Neither the window 
width I , nor replacing the "empty" window by an "opaque" one with a suitable 
point interaction as in Ref. 23 allow for a sensible perturbation theory, because in 
both cases the unperturbed bound state disappears as the coupling is switched out. 
Hence we rely again on a numerical analysis based on the mode-matching technique; 
the results confirm our expectation about the resonance character of the scattering 
and its dependence on parameters of the problem. 

There is one more interesting aspect. The corresponding classical system of 
coupled ducts is pseudo-integrable, its phase space being of genus three. Other 
systems of that type have been recently studied 24 ' 25 ; it was shown that their quantum 
couterparts exhibit a chaotic behaviour. One asks naturally whether a similar effect 
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Figure 1: Laterally coupled quantum waveguides 

can be observed here. To find the level-spacing distribution of the bound states, 
a very wide window is needed to produce a large number of eigenvalues. At the 
same time, the spectrum has to be unfolded, i.e., rescaled so that the mean spacing 
does not change along it. The result suggests spacing distribution is then sharply 
localized around a fixed value, hence there is no chaos. This is not surprising, since 
all the bound-state wavefunctions have transversally the shape of the first mode, co 
effectively they correspond to a one-dimensional system. What is less trivial is that 
the spacing distribution of the scattering phase shifts also does not witness of a fully 
developed chaos; this suggests that repeated reflections are an essential ingredient 
of the chaotic behavior of particles in bounded pseudo-integrable billiards. 

The above mentioned scattering analysis relies on the stationary approach. The 
time evolution of wave packets would deserve a separate study. In this paper we 
limit ourselves to a single example: in the concluding section we present a numerical 
method to solve the corresponding time-dependent Schrodinger equation, which 
allows us to draw some qualitative conclusions about time delay in the scattering 
on the connecting window. 

II Preliminaries 

The system we are going to study is sketched on Fig. 1. We consider a Schrodinger 
particle whose motion is confined to a pair of parallel strips of widths d±, d 2 , respec- 
tively. For definiteness we assume that they are placed to both sides of the x-axis, 
and they are separated by the Dirichlet boundary everywhere except in the inter- 
val (—a, a) ; we shall denote this configuration space by Q and £ := 2a . Putting 
h 2 /2m = 1 , we may identify the particle Hamiltonian with the Dirichlet Laplacian, 

H = H(d 1 ,d 2 ;£) := -A« , (2.1) 

on L 2 (f2) defined in the standard way — see Ref. 12, Sec. XIII. 15 — since the 
boundary of Q has the segment property, it acts as the usual Laplace operator with 
the Dirichlet condition at the boundary. 

A simple bracketing argument 12 shows that H has bound states for all £ large 
enough. Let us first introduce some notation. Set d := max{di,d 2 } and D := 
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d\ + d 2 , and furthermore 



min{(ii, d 2 } 



max{di, d 2 } 

We shall also use fid := (^j , and /i£ corresponding in the same way to D and 
£ , respectively. Cutting now Q by the additional Neumann or Dirichlet boundaries 
parallel to the y-axis at x = ±a , we get H { t N) © ifW < # < # t (D) © # C (D) , where 
the "tail" part corresponds to the four halfstrips and the rest to the central part 
with the Neumann and Dirichlet condition on the vertical boundaries, respectively. 

Since a ess (H^) = [/^,oo), j = N,D, the same is by the minimax principle 
true for H , and possible isolated eigenvalues of H are squeezed between those of 
j = N,D . The Neumann estimate tells us that 

inf a(H(d u d 2 - £)) > fi D = /i d (l + v)~ 2 . (2.2) 

On the other hand, H has an eigenvalue below fxjj provided does, which is 

true if hd + fJ>e ^ j this shows that a sufficient condition for H(di, d 2 ; £) to have 
at least a bound state is that the length of the opening satisfies the inequality 

d(l + v) 



I > 



[2.3) 



1.155 ; it grows as fl becomes 



If v — 1 , the coefficient on the right side is 2/ 
more asymmetric. 

More generally, the number of eigenvalues of is N D := y (fi d — (j, D )/ ' ^ 

where [•] denotes the entire part (recall that since D < 2d , the first transversal 
excited state is already above /id), while the number of the "Neumann" eigenvalues 
is N N := 1 + N D ; this means that the number of bound states of H(d 1: d 2 ;£) 
satisfies the inequality 



y 



£ ^v{v + 2) 
d l + i/ 



< iV < 1 



£ tJv{v + 2) 
d l + i/ 



(2.4) 



We see that H(d±, d 2 ; I) has isolated eigenvalues, at least for £ large enough, despite 
the absence of (a nonzero-measure set of) closed classical trajectories mentioned in 
the introduction. In heuristic terms, this may be understood as a manifestation 
of the fact that the semi-infinite "spikes" of the opened barrier between the two 
ducts are capable of reflecting a quantum particle due to a finite smearing of the 
wavepacket. In the same way, one finds that the mth eigenvalue ji m of H(di, d 2 ; £) 
is estimated by 

/™ 1 x 2 ,, -i ,_ N 2 

(2.5) 



m 



1 



< 



1 



< (- 



A J [id (1 - 

where A := £/d, and that the critical value \ m = £ m /d at which mth eigenvalue 
appears satisfies the bounds 



(m - 1)(1 + v) 
u(u + 2) 



<\ m < 



mil + is) 



1/(1/ + 2) 



(2.6) 
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To learn more about the dependence of the eigenvalues and the corresponding eigen- 
functions on A and v , we have to use a different technique. 



Ill Existence of bound states 



The above existence argument giving ( |2.3[ ) is a crude one; in fact, there is no lower 
bound on the window width as the following result shows: 

Theorem: H(di, d 2 ; t) has an isolated eigenvalue in \p, D , fi d ) for any £ > . 

Proof: We modify for the present purpose the variational argument of Ref. 6; see 
also Ref. 8, Sec. 2. Without loss of generality we may assume that d 2 < d\ = d. 
The transverse ground-state wavefunction is then 



xi(y) 



i sin ( K iy) ••• ye (MO 

... otherwise 

where K\ := yfjld ; similarly we define the transverse ground state in the opening, 

Vi(y) ■= ]J-p sin(iTi(d 1 -y)) 

with K\ := y/fj,D ■ For any $ G D(H) we put 

:= \\H®\\ 2 - fi d \\<f>\\ 2 

(if not marked explicitly, the norms always refer to L 2 (Q) ). 

Since the essential spectrum of H starts at fid , we have to find a trial function 
$ such that g [$] < ; it has to belong to the form domain Q(H) which means, in 
particular, that it must be continuous inside Q but not necessarily smooth. Notice 
first that if $(x,y) = (p(x)xi(y) , we have 

To make the longitudinal contribution to the kinetic energy small, we use an external 
scaling. We choose an interval J := [—b, b] for a positive b > a and a function 
ip G S(M) such that <p(x) — 1 if x G J ; then we define the family {{p a : a > } 
by 

!<f(x) ... \x\ < b 

<p(±b + a(x =F b)) ... \x\>b 
Finally, let us choose a localization function j e C^°((—a,a)) and define 

KA x ,y) : = <Mz)[xi(y) + £j(z)V(y)] (3.2) 
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for any a, e > . The main point of the construction is that we modify the factorized 
function we started with in two mutually disjoint regions, outside and inside the 
rectangle J x (— cf 2 , d\) . Hence the functions ip' a and j 2 have disjoint supports. 
Using this together with the identity 

II ' II 2 _ || '||2 

WAl^{]R) — a \\ L P \\L 2 (M) 

and the explicit forms of the functions Xh Vi > we substitute (|3.2|) into (|3.1|) and 
find after a tedious but straightforward computation 

q [^, £ } = awwhirn - 4^r 3/2 £~ 1/2 iij 2 ii! 2( iR) sin (rr^) 

(3.3) 

By construction, the last two terms on the right side of (|3.3| ) are independent of a . 
Moreover, the term linear in e is negative (recall that v £ (0, 1] ), so choosing e 
sufficiently small, we can make it dominate over the quadratic one. Finally, we fix 
this £ and choose a small enough a to make the right side of (|3.3j ) negative. | 

Remark: Though it is not the subject of the present paper, we want to note that 
the same argument demonstrates existence of a bound state in a straight Dirichlet 
strip with an arbitrarily small protrusion; one has only to replace J x [— d 2 , 0] by 
a rectangle contained in the protruded part. An alternative proof of this result has 
been given recently in Ref. 26; these autors also derived an asymptotic formula for 
the eigenvalue in terms of the protrusion volume. 



IV Mode matching 

To learn more about the eigenvalues and eigenfunctions in question, we shall now 
solve the corresponding Schrodinger equation numerically. Since Q consists of sev- 
eral rectangular regions, the easiest way to do that is by the mode-matching method. 



IV. 1 The symmetric case 



Consider first the situation when d\ = d 2 = d . The Hamiltonian (|2.1| ) then de- 
couples into an orthogonal sum of the even and the odd part, the spectrum of the 
latter being clearly trivial, i.e., the same as in the case I = . At the same time, 
the mirror symmetry with respect to the y-axis allows us to consider separately the 
symmetric and antisymmetric solutions. 

We may therefore restrict ourselves to the part of Q in the first quadrant, with 
the Neumann boundary condition in the segment (0, a) of the rr-axis, and Neumann 
or Dirichlet condition in the segment (0, d) of the y-axis. We expand the sought 
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solutions in terms of corresponding transverse eigenfunctions 

Xj(y) ■= y ^ sin^y) , J = 1,2,... (4.1) 

faiy) := V2 V2j ^(y) = ^| sin^^d - y)) , j = 1,2,... , (4.2) 

where Kj := j/ti and K^j-x '■= (2j — V)K\ . A natural Ansatz for the solution of an 
energy efi d , \ < t < 1 , is 

00 

1>{z,y) =Y,h^ a ~ x) xM (4-3) 
i=i 

for x > a , where ^ := Kiy/ j 2 — e , and 

~ cosh(p,x) , . . . , , ^ sinh(p,x) 
^ cosh(pja) pi smh(p i a) 

for < x < a and the symmetric and antisymmetric cases, respectively, where the 



longitudinal momentum is defined by pj := KiJ (j ' — ~) — e . It is straightforward 

to compute the norms of the functions ( |4.3|) and (|4.4| ); since j" 1 ^ and tend 
to Hd as j — > 00 , the square integrability of ip requires the sequences {a,-} and 
{bj} to belong to the space ^ 2 (j _1 ) • 

As an element of the domain of H, the function ip should be continuous to- 
gether with its normal derivative at the segment dividing the two regions, x — a. 
Let us first solve this condition formally. The continuity means Y^T=i a k < f ) k{y) = 
Hk^i^kXkiy) ; using the orthonormality of {xj} we get from here 

00 

bj = a k(Xj,<Pk) • (4.5) 

k=l 

In the same way, the normal-derivative continuity at x = a yields 

00 

<?A + H a k p k tanh(p k a) (xjAk) = (4.6) 
fc=i 

in the Neumann case, and the analogous relation with tanh replaced by coth for 
Dirichlet. Substituting from ( |4.5|) to (4.6 ), we can write the equation as 

Ca = 0, (4.7) 

where 

C jk ■= Uj+Pkl c ^ h \(Pka)) (XjAk) (4-8) 
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in the Neumann and Dirichlet case, respectively, with the two orthonormal bases 
related by 

= — (4.9) 

n f - (k - 1 



One has to make sure, of course, that the equation (|4.8|) makes sense, and 
that one can solve it by a sequence of truncations. It is possible to follow the 
procedure formulated in Ref. 4. A more direct way, however, is to notice that if ip 
is an eigenvector of H , it must belong to the domain of any integer power of this 
operator. It is easy to check that if) G D(H n ) iff {%}, {bj} G £ 2 (j 2n ~ 1 ) ; hence 
the sought sequences should belong to £ 2 (j s ) for all s > —1 . This fact also justifies 
a posteriori the interchange of summation and differentiation we have made in the 
matching procedure. 

Consider now the diagonal operator S r on t 2 (j~ l ) , (S r a)j := j~ r a,j . If C has 
zero eigenvalue with a fast decaying eigenvector, the same is true for C^'^ := S s CS r 
with arbitrary non-negative s, r . The last named operator is represented by the 
matrix _ ^ 

Cjk r) ■= (qj + Pk tanh(p fc a) ) ^— ^ 2j - 2 

71 J 2 -(k-k 



(for the sake of brevity, we speak about the Neumann case only) , so it is Hilbert- 
Schmidt for r, s large enough, and its eigenvalues can therefore be obtained from a 
sequence of truncated operators. Since finite matrices pose no convergence problems, 
the truncation procedure may be applied to the operator C directly. 

Of course, may have eigenvectors to which no square-summable eigenvec- 

tor of C corresponds, because S' 1 is unbounded for r > . Fortunately, the search 
for solutions may be terminated once we find the number of them which saturates 
the upper bound derived in Section 2. 

IV. 2 An alternative method 

A natural modification of the above described procedure is to express {a^} from 
(|4.5|) using the orthonormality of {(pk} , and to substitute it into ( (4.6[) ; then the 
spectral condition acquires the form 

b + Kb = 0, (4.10) 

where 

oo 

K jm := — (Xj,<Pk)Pktanh(p k a) (<frk,Xm) , (4.11) 

and the same with coth(pfca) in the Dirichlet case. 

The two approaches are, of course, equivalent. Solving the equation numeri- 
cally, however, we truncate not only the matrices but also the series in (|4 .11 ). The 



sequences approximating a given eigenvalue are therefore different. Moreover, in 
the examples given below we find them monotonous in the opposite sense. The 
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sequences coming from ( f4.7|) were approaching the limiting values from above, while 
those obtained from ( |4 . 1 0| ) were increasing; in combination this gives a good idea 
about the numerical stability of the solution. 



IV. 3 The asymmetric case 

Let us pass now to the case, when the widths of the ducts are nonequal, d\ ^ d,2 . 
Without loss of generality, we may again suppose that d 2 < d\ = d . With the 
mirror symmetry with respect to the y-axis in mind, we shall consider the right- 
halfplane part of Q only with the Neumann and Dirichlet condition on the segment 
W := [— d2,di] of the y-axis. 

To expand the sought solution, we need again suitable transverse bases. In the 
"connecting part" , < x < a , we use 

rik{y) = ^| sin^K -y)), k — 1,2,... , (4.12) 
where Kj, := kK\ = A;/t 1 (l + z/) -1 . On the other hand, for the ducts we choose 

xf\y) ■= Msm( Kj y)i + (y), j = l,2,..., (4.13) 

X [ f\y) ■= ~ \H- sinO^- 1 ^-^) , j = 1,2,..., (4.14) 

V "2 

where Kj := ji^i and i± are the indicator functions of the intervals T> + := [0,di] 
and D- := [—^2,0] , respectively. 

The union of the two bases is, of course, an orthonormal basis in L 2 (W) . Since 
the numerical computation involves a truncation procedure, we need to introduce 
a proper ordering. For that we arrange the eigenvalues corresponding to fl4.13|) 
to a single nondecreasing sequence. Equivalently, we arrange the numbers j, kv~ x 
with j, k = 1,2,... into a nondecreasing sequence (if v is rational and there is a 
coincidence, any order can be chosen in the pair); we denote its elements by 6 m , 

9 x :=l, 6 2 := min{2,z/~ 1 }, etc. 

The corresponding ordered basis in L 2 (W) is 

f x f\y) ••• o m =j 
U : U(y) = (4.15) 
1 x { f\y) ••• e m = jv- 1 

Consider first the even solutions, i.e., the Neumann condition at x = . A natural 
Ansatz for a solution of an energy e/i^ , (1 + v)~ 2 < e < 1 , is 

~ cosh(p fc x) 
£^ cosh(p fc a) 
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(4.16) 



^(z.y) := £&f } e^ (a " X) xf } (2/) ••• x>a;yGP ± , 
j'=i 



where 



and 



:= «a yjf - e , g< } := Kl ^ (£) - e 
The duct part of Q4.16 ) can be also written in a unified way as 



ij>(x,y) = £ c m e r -^U(y), (4-17) 



D r m (a-x) 
m=l 

where 

bf ] ... 0,„ j [ ... 0,„ j 

and r m : = < 

^ ••• 9 m =jv- 1 { ... B m = jv~ 1 

Using the continuity of the function and its normal derivative at x = a together 
with the orthonormality of {x^} > we find conditions for the coefficient sequences, 

oo 



bf = EMxfW), (4-18) 

k=l 

oo 

€ )h t ) + E a kPk tanh(p fc a) (xf\%) = 0. (4.19) 



k=i 

This can be also written as 



a fc p fc tanh(p fe a) (£ m ,rj k ) = ; 

fc=i fc=i 
substituting from the first equation to the second one, we obtain the spectral con- 



dition in the form (fO) with 

C mk := (r m + pfctanh (p fc a)) (£m,Vk) , (4.20) 



where the overlap integrals are given by 



(xS +) ,%) 



(Xj \vk) 



2j sin 



TTy/T+U j2 _ ( 



1+u 



2 ■ 



2 j TIT sin(^) 



(4.21) 



7T V 1 + V A2 _ ( JSK_ 

ll 



2 ' 



In the odd case, i.e., Dirichlet condition at x = , we get the same equation with 
tanh replaced by coth in ( |4.20| ). 

By a straightforward modification of the above argument, one can check that the 
coefficient sequences have a faster-than-powerlike decay and the spectral condition 
can be solved by a sequence of truncations. One can also rewrite the condition in 
the form analogous to ( |4.10| ) , c + Kc = , where 

J oo 

K jm ■= — Uj,Vk)Pk^Bh(p k a) (%,£ m ) . (4.22) 
r i k=i 

V Scattering 



The analysis is similar to that of the previous section. The incident wave is supposed 
to be of the form Xj (y) e~ lkj x in the upper channel, where we have introduced 



kf ] := k x Jk 2 -j 2 , fcH := ki \lk 2 - ( 3 - 



v 



we denote by rj^, £^7 , respectively, the corresponding reflection and transmission 
amplitudes to the j'-th transverse mode in the upper/lower guide. Due to the 
mirror symmetry, we can again separate the symmetric and antisymmetric situation 
with respect x = and to write 

r (±) _ I (M±) , («,±n A±) _ I ( («,±) _ («,±n t h 
' jj' ~ 2 V p ii' p 5i' ) ' jj' ~ 2 \ Pjj ' Pjj / ' 

where pff^ , a = s,a , are the appropriate reflection amplitudes. In the symmetric 
case we have the following Ansatz for the solution 

~ cos(p e x) 

Vfoi/) := 2^ a e 7 cW) ••• 0<a;<a, 

^ cos(p^a) 

V>(*,v) := E(%'e-^ +) (— )+ p W e ^ ) ^) ... x > a; y £ T> + , 

i'=i v y 



(5.2) 



oo 



2/) : = 2 Pj? e%ki> iX a) Xj> \y) . . . x > a; y eV_ . 



'33' 

f=X 



The last two relations can be written also as 

oo 



m'=l 

where 



Pjj' ■ ■ ■ 9m — 3 , Qm> — f 

k 
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Matching the functions (|5.2| ) smoothly at x = a we arrive in the same way as above 
at the equation 

oo 

(ih+ Vm> tan(p m /a) ) (& 2ikg5 m £ , (5.3) 

m'=l 

where the index m corresponds to the incident wave and the overlap integrals are 
given again by ( |4.21|) ; in the antisymmetric case on has to replace tan by —cot . 
The reflection amplitudes are given then by 

oo 

Pml = Srni + a m> (&,Vm') ! (5-4) 

m'=l 

they determine the full S-matrix via ( |5.1|) . 

VI The results 
VI. 1 Bound states 

The results of the mode-matching computation are illustrated on Figures 2-4. In 
accordance with the general results of Section 2 the eigenvalues decrease mono- 
tonously with the increasing window width and one can sandwich them between 
the estimates ( |2.5|) . The eigenfunctions decay exponentially out of the "interaction" 
region. The ground state wavefunction is, of course, positive up to a phase factor; 
the nodal lines of the excited states are parallel to the y-axis. The last feature 
illustrates once more that apart of the exponential tails in the ducts, the quantum 
particle "feels" the window part as a closed rectangular resonator. 

It is also interesting to estimate the rate at which the eigenvalues emerge from 
the continuum. The results of the mentioned paper 26 together with the Dirichlet 
bracketing allow us to find a simple upper bound for the ground-state energy by 
means of a single strip with a "blister" whose volume is squeezed to zero. Since 
the asymptotic formula derived in Ref. 26 applies to "gentle" protrusions, it may 
be employed if the power with which the bump is scaled transversally is larger than 
the longitudinal one. Hence the gap between the eigenvalue and the continuum for 
a narrow window is bound from below by C(e)£ A+£ + 0(£ 5 ) for any e > . 

This can be compared with the numerical results. Redrawing the first eigenvalue 
curve of Fig. 2 and analogous results for v ^ 1 in the logarithmic scale, we find 
that the asymptotic behavior is powerlike. The convergence of our method for small 
i is rather slow; nevertheless, using cut-off dimensions of order 10 3 we get for the 
power values witnessing clearly that the above bound is saturated, 

= fid- c{v)t + 0{t) . (6.1) 

The numerically found coefficient civ) is monotonous and reaches its maximum 
value for v — 1 ; this is the expected behavior as can be seen from a simple bracketing 
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argument. Proving the conjecture (|6.1| ) and finding an analytical expression for 
c{v) remains an open problem; the same can be said about the "coupling-constant 
thresholds", i.e., the way the other eigenvalues emerge from the continuum. 

VI. 2 Scattering 

The passage of the particle through the window region is determined by the trans- 
mission and reflection amplitudes ( |5.1|) . The physically interesting quantity is the 
conductivity. If we suppose, for instance, that the particle comes from the upper 
right guide and leaves through the upper left one, then the conductivity (denoted 
conventionally as TP and measured in the standard units 2e 2 /h) is given by 

Ik] L.C+) 

G(k) = £ T&yl^COr. (6.2) 

j,3'=l 

and similarly for the other combinations; the summation runs over all open channels. 
The resonance structure is visible on Figure 5. 

Another insight can be obtained by investigating the probability flow distribu- 
tion associated with the generalized eigenvector ( |5.2|) which is defined in the standard 
way, 

j(x) := —ii/}(x)Vtp(x) . (6.3) 

The flow patterns change with the momentum of the incident particle. They ex- 
hibit conspicuous vortices at the resonance energies which represent the "trapped 
part" of the wavefunction; this phenomenon is illustrated on Figure 7. It has been 
argued in the literature that leaky wires similar to those studied here may serve 
as switching devices 22 . The vortices which emerge in resonance situations lead to 
the appearance of a magnetic dipole moment, which might be in principle measured 
experimentally. In this respect situations with a single well developed vortex such 
as the one illustrated on Figure 6 are particularly promising. 

VI.3 Chaos 

Discussing a chaotic behavior of a quantum system, it is useful to start with its 
classical counterpart, and in particular, its phase space. In the present case of an 
infinite two-strip "billiard" there are no closed classical trajectories with exception of 
the obvious zero-measure set, hence one has to consider the scattering, i.e., motion 
of a point particle bouncing its way through the system; the reflection from the walls 
is supposed be perfectly elastic. 

There are two integrals of motion: the longitudinal component of the momen- 
tum, Ii = p x , and the modulus of its transverse part, I<i = \p y \ ■ Hence the phase 
space trajectory of the system is restricted to a two-dimensional manifold (invari- 
ant surface) in the four-dimensional phase space. However, due to the singularity 
of corresponding classical flow at the edges of the connecting window, the topology 
of this surface is not equivalent to that of a two-dimensional cylinder, but rather 
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of a pair of mutually crossed cylinders; similar systems are usually dubbed pseudo- 
integrable 27 . The topological structure of the invariant surface has a consequence 
for the quantum counterpart: the system cannot be quantized semiclassically. 

On the other hand, the quantum system of coupled waveguides has in view of our 
previous arguments bound states, even many of them iff i^fv ^> d . Then one can 
plot the distribution of the eigenvalue spacing as shown on Figure 8 for a particular 
value of v ; the character of the distribution does not change as the latter is varied. 
The natural unfolding means in this case to employ the corresponding momentum 
value ipi = (e — (l + z/)~ 2 ) . The results differ from typical (unfolded) eigen- 
value distributions in billiards, both integrable and chaotic, in the first place due 
to the existence of the sharp localization around a fixed value. The used statistics 
(several thousand eigenvalues) does not allow us to tell what is the behavior around 
zero; we see, however, that the decay off the peak is at least exponential. This differs 
substantially from a typical behavior of chaotic systems, however, one should not 
be surprised because all the corresponding eigenfunctions are dominated transver- 
sally by the lowest mode, so the bound-state family in our "billiard" is effectively 
one-dimensional . 

It is less trivial whether a chaotic behavior may be manifested in the scattering; 
recall that spatially restricted pseudo-integrable billiards are known to exhibit the 
so-called wave chaos 24 . To decide whether a quantum scattering system is chaotic or 
not, one has to study eigenvalue distribution of the corresponding S-matrix, again 
properly unfolded, which is expected to conform with that of the Dyson circular 
ensemble of random matrices 28 in the former case. We have performed this task 
for the system under consideration numerically, analyzing the distribution of the 
spacing between two neighboring eigenvalues of the S-matrix. The result is plotted 
on Figure 9; they are compared with the Wigner and the Poissonian distributions 
peculiar for the chaotic and non-chaotic situation, respectively. It can be seen that 
the overall shape of this distribution matches the Poissonian distribution for all 
spacings large enough; on the other hand, the deformation of the distribution near 
the origin provides a clear sign of non-integrability of the system. The fact that 
this non-integrability differs from a typical chaotic behavior can be attributed to 
the fact that the scattered particle passes the window region "only once" without 
being bounced to and fro as it is the case of finite billiards. 

The absence of the fully developed chaos in the coupled waveguides can also 
be seen when plotting the coefficients an which determine the wavefunction in the 
interaction region by (|5.2| ) as illustrated on Figure 10. Their distribution remains 
well localized even for higher energies of the incoming particle, its tail being ap- 
proximately exponential, while in case of an irregular scattering one would expect a 
slower decay. 
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VII Time evolution 



Up to now we have discussed the coupled waveguide system from the stationary point 
of view only. Let us look briefly how the window coupling can affect propagation 
of wavepackets in the ducts. This problem has a natural motivation: it has been 
suggested recently 22 ' 29 ' 30 that coupled electron waveguides provide an analogue of 
the optical directional coupler in the sense that they may switch electrons from one 
quantum wire to another. Moreover, the authors of Ref. 22 conjectured that the 
electron switching process should be rather fast due to the direct character of the 
corresponding resonance, since the electron is not trapped in the interaction region 
during the resonant switching. 

The existence of probability-flow vortices discussed above in the interaction 
region indicates that this might not be the case, i.e., that the electron dwelling time 
in the junction may not be generally neglected. To get a better insight we have 
investigated time evolution of wave packets numerically. This can be achieved by 
approximating the evolution operator by a Trotter-formula product — see Ref. 12, 
Sec. VIII. 8 — with the Dirichlet boundary condition replaced by a very steep and 
narrow potential barrier localized along the boundary; the latter has been chosen in 
such a way that the dynamics of the system was equivalent to the dynamics of the 
true Dirichlet problem for all times taken into account, i.e., that the tunneling leak 
was negligible during that period. 

The kinetic- and potential-part factors of the evolution operator are then mul- 
tiplication operators in the momentum and coordinate representation, respectively; 
the passage between the two representations has been realized by means of the two- 
dimensional Fast Fourier Transform method 31 with a grid of 2 9 x 2 7 points. The time 
evolution of a wavepacket approaching the junction through the upper right arm of 
the structure is plotted on Figure 11. The incoming wavefunction was chosen as 
ip(x,y) := g(x)xi (y) , where g(x) := exp{— a(x— x ) 2 + ikx} with suitably chosen 
parameters a, xo ■ 

The difference between the resonant and nonresonant situation is clearly visible. 
In the first case the electron stays in the junction region and escapes only slowly, 
while the electron whose momentum is localized around a slightly different but non- 
resonant value of momentum passes the junction "ballistically" . Wang and Guo 22 
based the mentioned conjecture — which in a realistic situation would lead to ul- 
trashort switching times of a few ps only — on a concept of transmitivity of coupled 
waveguides leaning on a classical intuition. As we have said in the introduction 
and demonstrated in the previous sections, this may be a false guide when quantum 
systems are considered. The example of time evolution offers another illustration. 
During the resonance-scattering process the evanescent-mode amplitudes inside the 
quantum wire are considerably enhanced; as a result the electron is trapped tem- 
porarily inside the juction. The probability of finding it there in the resonant and 
nonresonant case, respectively, is shown on Figure 12. It is desirable to perform 
the time-delay analysis for the present model, in particular, to confirm that the 
"switching time" of the coupler is inversely proportional to the resonance width. 



16 



Acknowledgement 

The work has been partially supported by the Grants AS No. 148409 and GA CR 
No.202-93-1314. 

References 

[1] R. Hempel, L.A. Seco, B. Simon: "The essential spectrum of Neumann 
Laplacians on some bounded singular domains", J. Funct. Anal. 102, 448- 
483 (1991) 

[2] B. Simon: "The Neumann Laplacian of a jellyroll", Proc. Amer. Math. Soc. 
114, 783-785 (1992) 

[3] P. Exner, P. Seba: "Bound states in curved quantum waveguides" , J. Math. 
Phys. 30, 2574-2580 (1989) 

[4] P. Exner, P. Seba, P. Stovicek: "On existence of a bound state in an 
L-shaped waveguide", Czech. J. Phys. B39, 1181-1191 (1989) 

[5] R.L. Schult, D.G. Ravenhall, H.W. Wyld: "Quantum bound state in classi- 
cally unbound system of crossed wires", Phys. Rev. B39, 5476-5479 (1989) 

[6] J. Goldstone, R.L. Jaffe: "Bound states in twisting tubes" , Phys. Rev. B45, 
14100-14107 (1992) 

[7] G. Dunne, R.L. Jaffe: "Bound states in twisted Aharonov-Bohm tubes", 

Ann. Phys. 233, 180-196 (1993) 
[8] P. Duclos, P. Exner: "Curvature-induced bound states in quantum wave- 
guides in two and three dimensions", Rev. Math. Phys. 7, 73-102 (1995) 
[9] P. Exner: "A model of resonance scattering on curved quantum wires", 
Ann. Physik 47, 123-138 (1990) 
[10] P. Duclos, P. Exner, P. Stovicek: "Curvature-induced resonances in a two- 
dimensional Dirichlet tube", Ann. Inst. H. Poincare 62, 81-101 (1995) 
[11] P. Duclos, P. Exner, B. Meller: "Exponential bounds on curvature-induced 
resonanaces in a two-dimensional Dirichlet tube" , submitted for publication 
[12] M. Reed, B. Simon: Methods of Modern Mathematical Physics, I. Func- 
tional Analysis, IV. Analysis of Operators, Academic Press, New York 1972, 
1978. 

[13] CM. Savage, S. Markensteiner, P. Zoller: "Atomic waveguides and cavi- 
ties from hollow optical fibres", in Fundamentals of Quantum Optics III 
(S. Eklotzky, ed.), Springer, Berlin 1993. 

[14] P. Exner, P. Seba: "Trapping modes in a curved electromagnetic waveguide 
with perfectly conducting walls", Phys. Lett. A144, 347-350 (1990) 

[15] J. P. Carini, J.T. Londergan, K. Mullen, D.P. Murdock: "Bound states and 
resonances in quantum wires", Phys. Rev. B46, 15538-15541 (1992) 

[16] J. P. Carini, J.T. Londergan, K. Mullen, D.P. Murdock: "Multiple bound 
states in sharply bent waveguides", Phys. Rev. B48, 4503-4514 (1993) 

[17] B. Simon: "The bound state of weakly coupled Schrodinger operators in 
one and two dimensions", Ann. Phys. 97, 279-288 (1976) 



17 



[18] R. Blanckenbecler, M.L. Goldberger, B. Simon: "The bound states of 

weakly coupled long-range one-dimensional quantum Hamiltonians" , Ann. 

Phys. 108, 69-78 (1977) 
[19] Y. Hirayama, Y. Tokura, A.D. Wieck, S. Koch, R.J. Haug, K. von Klitzing, 

K. Ploog: "Transport characteristics of a window-coupled in-plane-gated 

wire system", Phys. Rev. B48, 7991-7998 (1993) 
[20] Ch. Kunze: "Leaky and mutually coupled quantum wires", Phys. Rev. 

B48, 14338-14346 (1993) 
[21] LA. Larkin, E.V. Sukhorukov: "Method to investigate the random potential 

in a quantum point contact", Phys. Rev. B49, 5498-5507 (1994) 
[22] Jian Wang, Hong Guo: "Time-dependent transport in two-dimensional 

quantum-wire structures", Phys. Rev. B48, 12072-12075 (1993) 
[23] P. Exner, M. Tater: "A one-band model for a weakly coupled quantum- 
wire resonator", Phys. Rev. B50, 18350-18354 (1994) 
[24] A. Shudo, Y. Shimizu, P. Seba, J. Stein, H.J. Stockmann, K. Zyczkowski: 

"Statistical properties of spectra of pseudointegrable systems" , Phys. Rev. 

E49, 3748-3756 (1994) 
[25] K. Zyczkowski: "Classical and quantum billiards", Acta Phys. Polonica 23, 

245-270 (1992) 

[26] W. Bulla, F. Gesztesy, W. Renger, B. Simon: "Weakly coupled bound 
states in quantum waveguides", submitted for publication 

[27] P.J. Richens, M.V. Berry: "Pseudointegrable systems in classical and quan- 
tum mechanics", Physica D2, 495-512 (1981) 

[28] R. Blumel, U. Smilansky: "Classical irregular scattering and its quantum 
mechanical implications", Phys. Rev. Lett. 60, 477-480 (1988) 

[29] J. A. del Alamo, C.C. Eugster: "Quantum field-effect directional coupler", 
Appl. Phys. Lett. 56, 78-80 (1990) 

[30] N. Tsukada, A.D. Wieck, K. Ploog: "Proposal of novel electron wave cou- 
pler devices", Appl. Phys. Lett. 56, 2527-2529 (1990) 

[31] D. Kosloff, R. Kosloff: "Absorbing boundaries for wave propagation prob- 
lems", J. Comput. Phys. 63, 363-376 (1986) 



Figure captions 

Figure 1 Laterally coupled quantum waveguides 

Figure 2 Bound-state energies vs. the window width £ in the symmetric case. 

Figure 3 The ground-state eigenfunction in the symmetric case for £/d — 0.3 . 

Figure 4 The eigenfunction of the second excited state in the unsymmetric case, 
v = 1/2 Jot £/d 1 = 1.08. 
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Figure 5 The conductivity for the particle coming from the right in the upper 
duct as a function of the momentum k and the width g?2 of the lower tube 
for d± = 7T, £ = 2 . (a) The particle leaves through the upper left channel. 
A deep resonance is clearly visible, (b) The particle leaves through the lower 
left channel. The conductivity is zero when there are no propagating modes 
in the lower part. 

Figure 6 The quantum probability flow ( |6.3| ) for the symmetric situation, v = 
1 , in the resonance and non-resonance situation, respectively. The appearance 
of vortices associated with the resonance scattering is obvious. 

Figure 7 A single vortex corresponding to the sharp stopping resonance of 
Figure 6a. The conductivity is small in this situation so the waveguide system 
is closed for the electron transport. 

Figure 8 The unfolded level-spacing distibution of the symmetric and antisym- 
metric bound states for v = 2(1 + Vb)^ 1 . 

Figure 9 The unfolded level-spacing distribution for the S-matrix correspond- 
ing to v = 2 /it and averaged over momentum, in comparison with the Poisson 
and Wigner distribution. 

Figure 10 The absolute value of the coefficients a>i of eq.(5.2) in the symmetric 
case for v = 2/tt and k = 28.432 ; the particle is supposed to be initially in 
the 18-th transverse mode. 

Figure 11 The time evolution of the wavepacket inside the junction with v = 2/n 
plotted for times t — 0, 5, 10, 15, and 20, respectively, (a) The resonance case 
with k = 1.4242, (b) the near-to-resonance situation, k = 1.48. 

Figure 12 The probability that the electron will be found within the junction as 
a function of time evaluated for the same parameters as on Figure 11. 
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